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CAS A X-RAY SPECTRUM: 

EVIDENCE FOR IRON LINE EMISSION 

P. J. Serlemitsos, E. A. Boldt, S. S. Holt, R. Ramaty 
NASA, Goddard Space Flight Center 

A. F. Brisken 
University of Maryland 


ABSTRACT 

A sensitive measurement by rocket borne detectors of the x-ray 

flux from Cas A has revealed a steep continuum and a broad spectral 

feature in the region where line radiation from iron nuclei would be 

-2 -1 

expected. The flux in this feature is .0122 + .0017 photons cm s ; 

-9 -2 -1 

the total energy flux from 2 to 10 keV is 1.02 x 10 ergs cm s . 

The presence of broad iron lines is consistent with a model in which 
~ 13 MeV/nucleon iron nuclei charge exchange with surrounding inter- 
stellar oxygen and other heavy atoms. The model suggests that a 
substantial fraction of the energy from the outburst has gone into low 
energy cosmic rays which produce the observed HII region surrounding 


the remnant . 
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I. INTRODUCTION 

Cas A is believed to be the youngest known supernova remnant in 
the galaxy resulting from an explosion that took place around 1667 A.D., 
about 2.8 kpc away (Van Den Bergh & Dodd 1970). Optical studies of 
the remnant reveal a rapidly expanding nebula surrounded by a faint 
HI I region which extends out to 5.7 pc from the site of the explosion. 
The remnant is a strong source of nonthermal radio emission being by 
far the brightest radio source in the sky. 

Byram, Chubb and Friedman (1966) first detected Cas A in the X-ray 
band in a 1965 rocket survey. Gorenstein, Kellogg and Gursky (1970) 
produced the first spectrum for this source in the band 2 to 10 keV. 

The UHURU catalogue of X-ray sources (Giacconi et al. 1972) lists 
Cas A as a steady emitter and one of the relatively few sources with 
which a definite optical identification has been made. The lack of 
temporal variations (Holt et al. 1973), as well as the nature of the 
optical and radio emissions strongly implies that the X-ray source is 
an extended object, much as the source in the Crab Nebula and unlike 
such thermal sources as Sco X-l and Cyg X-2, or the eclipsing binaries 
Her X-l and Cen X-3. 

In this paper we shall present new spectral information pertaining 
to the X-ray flux from Cas A, obtained from a rocket borne experiment 
launched from White Sands, New Mexico at 0900 UT on 19 May 1972. 

II. THE EXPERIMENT 

We scheduled most of the available time in the rocket flight to 
conduct a sensitive study of Cas A and of the nearby Tycho remnant. 
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Results from a temporal analysis of the data from both sources have 

been reported elsewhere (Holt et al. 1973). The data we are presenting 

here were obtained with a xenon-filled, multi-anode, multi-layer 

proportional counter, equipped with a .001" kapton window and a 3-degree 

o 

circular Cu collimator. The total open area is 711 cm . The exposure 
to Cas A lasted 61.4 seconds. Star field pictures showed that the 
detector was pointed to within 0.1 degrees of the source throughout 
the exposure, during which time the counting rate averaged 93.4 s -1 . 

The average background rate was measured at 16.8 s ^ . An onboard Fe'*'’ 
source provided calibration data before the doors to the experiment 
opened, as well as after they closed at the end of the flight. 

III. ANALYSIS 

We will present the data in a representation that, we believe, 
best allows for comparisons with other observations. We have computed 
an incident spectrum S(E) at energy E as follows: 


S(E) 


N obs (E) 

e(E) 


( 1 ) 


where N 0 ^ s are the net observed counts from the source, and e is an 
effective efficiency defined as 


e(E) 


N (E) 
calc 

S'(E) 


C f° R (E , E ’ ) S ' (E ' ) dE ' 
J o 

S ' (E) 


( 2 ) 


where R(E,E') is the detector response function as determined from 
source calibration, which relates the apparent energy to the incident 
true energy, S'(E) is a first approximation to the source spectrum, 
and C is a normalization constant. A usual procedure for arriving at 
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an approximate expression for S(E) has been to use equation 2 for 

2 

computing N ca i c (E)> which is then compared with N Q ^ S (E) using the y 
method. This may require a prohibitive number of trial functions 
when S(E) involves many parameters. In this paper we proceed 
to obtain better approximations for S(E) using the above integral 
equation, starting with S'(E) = constant. This procedure converges 
rapidly, resulting in a near final representation for S(E) even after 
only one iteration. It is, of course, obvious that any sharp features 
of the incident spectrum, such as lines, will be smoothed over by an 
amount commensurate with the energy resolution of the instrument. 

IV. RESULTS 

In Figure 1 we show the incident spectrum for Cas A inferred from 
an analysis as described above. No statistically significant net flux 
was detected above about 15 keV. The spectrum appears to be made up 
of more than just a single component. If we exclude the region 5 to 
10 keV, a good fit is obtained using a single power law with some 
absorption at the lower energies. The index to the best fit number 
spectrum is 4.5; the absorption implies interstellar matter in the line 
of sight of 1.4 x 10 hydrogen atoms cm , based upon the cross sections 
given by Brown and Gould (1970). We find no major segment of the data 
that can be similarly fitted' with a simple thermal continuum. Assuming 
that a power law with an additional component in the range 5 to 10 keV 
does adequately represent the incident spectrum, we have subtracted 
from the data the appropriate power law and we have plotted the residual 
flux in Figure 2. It is particularly significant that the broad 
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residual is centered near 7 keV, a region where the dominant K-shell 

+25 +24 

transitions from Fe and Fe would be expected. However, the 

feature in the data is too broad to be attributed to collisionally 
stripped ambient iron. Assuming that the residual in Figure 2 is due 
to iron lines, the observed broadening is not likely to have resulted 
from electron scattering of the radiation in a compact thermal source, 
but is probably due to electron cascades to the ground state in fast 
iron nuclei following charge exchange with neutral matter (cf. Silk and 
Steigman 1969). We note that the magnitude of the emission feature in 
relation to the continuum is prohibitively high (cf. Felten, Rees and 
Adams 1972) for a thermal spectrum. Furthermore, we reiterate that we 
could not fit any segment of the observed spectrum with a thermal con- 
tinuum, especially one with the high temperature needed to account for 
the presence of iron lines. Finally, we note that much of the previous 
evidence supports an extended source model for Cas A. 

If the interstellar region surrounding the remnant is permeated 
with low energy cosmic rays (i.e. a few MeV/nucleon) ejected by the 
supernova (cf. Ramaty et al. 1971) we can expect high Z nuclei to slowly 
lose their effective charge through accumulation of atomic electrons 
via charge exchange with ambient neutral atoms. A substantial number 
of transitions to the ground state from higher excited states will result 
in 2p- Is transitions at ~ 7.0 keV and 6.7 keV from H-type and He-type 
Fe respectively. These lines will subsequently be Doppler shifted to a 
degree depending on the velocity of the emitting ion. The final 
distribution will be a composite of all , contributing lines. Although 
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the contribution from lower Z elements could be comparable, the main 
effect will occur below about 2 keV where the total intensity is more 
than an order of magnitude greater than at 7 keV. 

To test this model, we fitted our data with the spectrum 


_-(E-6 .9) 2 /2ct 2 -j e -xCT g 


f = Ci[E- + C 27 ^ a e 


(3) 


22 -2 

where x = 1.4 x 10 hydrogen atoms cm and CJg is the effective absorp- 
tion coefficient for the interstellar medium (Brown and Gould 1970). 

2 

The best fit gives X = 17 for 12 degrees of freedom, and was obtained 
with 0/ = 4.5+.07; o = 1.23+.09 keV; Cj = 8.49 photons cm ^s ''"keV 

O 

C 2 = .00147. The errors shown are those that cause x to increase by 
1 when a corresponding parameter alone is allowed to change. If we 
set 


a CT intrinsic + 0 instrument (4) 

and substitute the measured detector resolution o^gtrument = .54 keV, 
we obtain CjL n trinsic = 1-1 keV. If we now assume that the 6.7 and 
7.0 keV lines contribute equally, we conclude that single lines at the 
source are broadened by about 33%. The corresponding particle energy 
is 13 MeV/nucleon. 

The residue in Figure 2 amounts to a flux of .0122+. 0017 photons 
-2 -1 

cm s . In the range 2 to 10 keV, we obtain a number flux of .223 

-2-1 -9 -2-1 

photons cm s and an energy flux of 1.02 x 10 ergs cm s . 

V. DISCUSSION 


Although Gorenstein, Kellogg and Gursky (1970) fitted the Cas A 
spectrum with a single power law they noted that the emission could be 
more complex since they did not obtain an acceptable fit. An examination 
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of their published data points indicates that the present model could 
provide a better fit. The amount of interstellar absorption inferred 
from their data is in good agreement with that given in this paper, 
both being consistent with 21-cm absorption measurements (cf. Davies 
and Shuter 1963). The UHURU catalogue gives an energy flux for Cas A 
of 9.1 x 10” ^ ergs cm’^s"^ in the range 2 to 10 keV, in excellent 
agreement with our result. 

A recent supernova remnant is a likely place to look for effects 
attributable to fast charged particles. If an appreciable fraction of 
the energy in the outburst goes into such particles, their effects may 
include the ionization of the surrounding region and the radiation which 
results from interactions with ambient matter. In ~ 300 years, the age 
of the Cas A remnant, 13 MeV/nucleon nuclei will traverse in straight 
line trajectories a distance of 15 pc, but the actual rectilinear distance 
away from the explosion may be smaller because of magnetic fields. 

Peinbert and Van Den Bergh (1971) discuss an HII region which surrounds 
the Cas A remnant and which extends out to 5.7 pc. From the emission 
measure they estimate an r.m.s. electron density in the region of 15 cm” . 
We shall explore the possibility of a closer tie between this observation 
and our charge exchange model. 

Nikolaev (1965) gives the following approximate expression for the 
cross section for charge exchange as a function of the velocity v of 
a highly stripped nucleus of charge Z incident upon a medium of atomic 


number Z me( j: 
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i 2 z 2(Z) 3 z 5 /V_) 8 [ (I_) 4 + 2 (l-) 2 ( Z 2 + id)+(Z 2 

o v r/ ™e d vo v Vc / v Vc / med n 2 med 


Z^ 2 i-5 


where a Q = .528 x 10" 8 cm, the radius of the first Bohr orbit, and 
v = 2.19 x 10 8 cm s _1 , the electron velocity in the first Bohr orbit; 
n is the principal quantum number of the state in which the electron 
is captured. In evaluating the cross section we summed over n > 1. 
Comparisons with experimental charge exchange cross sections (Heckman, 
Hubbard and Simon 1963; Nikolaev 1965; Macdonald and Martin 1971) 


indicate generally good agreement for v ^ Z mec j V o° ’'‘ n re ®'*' on 
~ 13 MeV/nucleon this condition is satisfied for all abundant inter- 


stellar species. In Table 1 we list cross sections calculated for 
various target elements. We also list the expected contribution to 
the total effective cross section per interstellar hydrogen atom based 
on universal abundances (Cameron 1967). Note that half of the total 


effective cross section is due to 0, and that only a small fraction 


may be attributed to light elements. This implies that charge exchange 
proceeds almost independently of the state of ionization of the medium 
up to temperatures of ~ 10 °K. 

The iron nuclei which produce the observed spectral feature could 


have reached charge equilibrium in the early stages of the expansion. 

At ~ 13 MeV/nucleon the equilibrium populations of Fe +28 , Fe +2 ^ and Fe +2 ^ 


are comparable (Heckman, Hubbard and Simon 1963) and hence iron line 
emission will proceed essentially as described above. At lower energy 
(< 10 MeV/nucleon) equilibrium charge states with bound K electrons 
are favored and, therefore, such line emission is suppressed. 
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A competing process involves the radiative capture by the fast 
nuclei of ambient or loosely bound electrons (cf. Schnopper et al. 

1972) which results in spectral features different from those we have 
discussed. An estimate, however, shows that the cross section for this 
process is more than an order of magnitude below that for charge 
exchange . 

Table 1. 

Charge Exchange Cross Sections Calculated for 13 MeV/nucleon 
Fe Nuclei Impinging on Interstellar Matter 


Element 

O 

a(cnr/ 

atom) 

— 2 

a (cm /hydrogen 

H 

1.2 

X 

io“ 22 

1.2 

X 

-22 

10 

He 

3.6 

X 

io" 21 

2.9 

X 

io' 22 

C 

5.7 

X 

io" 19 

2.9 

X 

IO" 22 

0 

1.7 

X 

10" 18 

15 

X 

10‘ 22 

Ne 

3.2 

X 

IO' 18 

2.9 

X 

io“ 22 

Mg 

4.8 

X 

IO" 18 

1.9 

X 

io- 22 

Si 

5.8 

X 

io' 18 

2.3 

X 

Hf 22 

Fe 

1.7 

X 

10" 18 

.7 

X 

io' 22 


Total a = 3 x IO -21 

The time for single electron capture by Fe nuclei may now be 
estimated from Table 1 as 2200/n^ years, where n^ is the total (neutral 
and ionized) interstellar hydrogen density. From our data we infer a 
rate of Fe line emission at the source of 1.1 x 10^ 8 photons s"'*'. 
Assuming one photon per Fe nucleus, the total number of iron nuclei is 

C O 

Np e = 7x10 /n^. For cosmic-ray abundances (Shapiro and Silberberg 



9 


1970), the corresponding number of hydrogen nuclei is = 1.8xl0^^/n^, 
so the total mass in fast particles is ~ 2 M^/n^, and the total energy 

C O 

is ~ 4x10 /n H ergs. If all the iron in the galaxy is similarly produced, 

Q 

we would require about 5x10 /n^ supernovae. 

Protons will be primarily responsible for the cosmic ray produced 
ionization in the HII region. The average electron density may be 
written as 


n„ 


M dE 

- N H x dt x T 
36xlO“ 6 V 


( 6 ) 


— - 1.46 x lO -12 !^ E - ' 3 MeV s" 1 


(7) 

-3 


where T = 10 10 s , the age of Cas A, and V = 2xl0 58 cm 3 , the volume of 
the HII region. The energy loss rate of a proton of energy E in neutral 

hydrogen is (Ramaty et al. 1971) 

dE 
dt 

Substituting (7) into (6) for E = 13 MeV we obtain n g — 17 cm This 

value is comparable to the r.m.s. density in the HII region. The 

coincidence, however, may be fortuitous because of a nonuniform density 

-3 

distribution. If the average density is less than 17 cm , most of the 
particle energy went into heat and expansion rather than ionization. 

The model thus predicts that essentially all of the interstellar gas 
in the nebula surrounding the remnant has been ionized by energetic 
particles which are now heating the gas. (However, charge exchange 
still takes place with the almost neutral heavier atoms.) A necessary 


consequence of this heating is bremsstrahlung x-radiation arising from 
collisions of suprathermal nuclei with ambient electrons (Boldt and 
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Serlemitsos 1969). If we assume that a cosmic ray abundance (Shapiro 
and Silberberg 1970) of energetic nuclei (relative to iron) were ejected 
at ^ 13 MeV/nucleon, then about 20% of the total x-ray flux at ^ 7 keV 
may be accounted for by suprathermal bremsstrahlung. Accounting for 
the complete continuum spectrum down to ~ 2 keV by this mechanism 
would require an order of magnitude more nuclei ejected at lower 
energies (~ 4-13 MeV/nucleon). 

The Crab Nebula is roughly 3 times as old as Cas A so that it too 
should be source of broad iron lines, if the two supernovae were 
similar events. The Crab x-ray spectrum has been extensively studied 
and, in the range of interest, it shows no evidence of spectral features 
(Boldt et al. 1969). We note, however, that the line flux that we have 
inferred for Cas A would at most amount to 2% over the strong and flat 
Crab continuum, and should not be easily detectable. Data from Tycho 
obtained in the same flight are now being analyzed. 

A test of our interpretation of the Cas A spectrum will be possible 
when more refined measurements of this source will be made. Specifically, 
the spectral feature must come from a spatially extended region and it 
should not show variability on time scales shorter than those of the 
supernova remnant. 
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FIGURE CAPTIONS 

Figure 1 - Incident x-ray spectrum from Cas A. The most accurate 

measure (with the exception of possible systematic errors) 

- of the differential photon flux is between 2.3 and 3.0 keV, 
and is .083+. 002 photons cm’V 1 keV _1 . 

Figure 2 - Residual flux after subtracting from the Cas A incident 

spectrum a power law spectrum with index 4.5 which includes 
interstellar absorption normalized to 1.4 x 10 22 hydrogen 
atoms cm . 
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